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Abstract 

We consider the effect of gluon radiation on the energy of top / antitop quarks and on 
the anticollinearity of top-antitop quark pairs produced in e + e~ annihilation. Our 
results are presented in terms of the ^-dependence of the tig cross section and the 
dependence on the cosine of the opening angle #12 between top and antitop for a 
center of mass energy of \f~cp' = 500 GeV . We then go on to determine mean values 
for the top quark's energy as well as its longitudinal and transverse projections, and 
for the deviation of sin 612 and cos #12 from the anticollinearity limits sin #12 = and 
cos #12 = —1. For a center of mass energy of 500 GeV we obtain (E q ) = 248.22 GeV , 
(E L ) = 247.24 GeV and (E T ) = 4.70 GeV. Thus, at this energy gluon radiation 
causes a total average energy loss of 0.71% of the top quark's energy. The average 
energy loss in the longitudinal direction is 1.06% and the average energy gain in 
the transverse direction is 1.88%. These percentage figures go up to 3.77%, 5.19% 
and 6.06%, respectively, at 1000 GeV. For the mean of the acollinearity angle O12 = 
180° - 6>i2 we obtain (6 12 ) = 1.25° at 500 GeV, the value of which goes up to 4.62° 
at 1000 GeV. From an analysis of the transverse momentum of the top we find that 
the mean transverse momentum of the top stays close to the mean total momentum 
of the gluon in the energy range from threshold to 1000 GeV showing that the gluon 
momentum has a large mean transverse component in this energy range. 
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1 Introduction 



Top quark pairs produced in e + e~ annihilation provide a very clean and kinematically 
well-defined source of top quarks. To leading order in the QCD coupling constant a s the 
energy of the top quark is exactly given by one-half of the center of mass energy. The 
top and antitop quarks decay so fast that there is no visible primary track. Even if the 
energy of the top quark is constrained by the beam energy one has to rely on kinematical 
reconstruction from an analysis of the decay products in order to pin down the momentum 
direction of the primary top quark pairs. Since to lowest order of QCD the top and antitop 
are produced exactly back-to-back, the lowest order anticollinearity provides for a further 
important kinematical constraint on the production kinematics. 

An exact knowledge of the momentum configuration of the primary produced top quark 
pairs is important for measurements of the spin observables of the top quark. If the top 
quark's momentum and momentum direction is known, one can boost to the rest frame of 
the top quark and do the spin analysis in the top quark rest system which is the optimal 
frame from the point of view of a spin analysis. As concerns spin-spin correlation effects, 
the anticollinearity of the produced top quark pairs is also important for optimizing spin- 
spin correlation measurements M. 

Previously it has been emphasized that the bremsstrahlung of gluons off the top quark 
pairs can have important effects on various characteristics of the it final state. For example, 
the mass determination of the top quark and the energy spectrum of secondary leptons 
will be affected by the presence of gluon bremsstrahlung R |3|, M. In order to quantify 
the effect of gluon bremsstrahlung, the authors of have calculated the gluon energy 
spectrum in e + e~ — > tig as well as the mean energy of the gluon (E g ). Having in mind 
the above-mentioned spin analysis we continue in the study of gluon emission effects in 
e + e~ — > tig and determine how the energy of the top (or antitop) quark is affected by 
gluon emission and how the anticollinearity of the top-antitop configuration is distorted 
by gluon emission. 

Our paper is structured as follows. In Sec. 2 we study the top energy spectrum and 
compare it to the bottom energy spectrum. From this we determine the mean deviation 
of the top and bottom energy E q from their lowest order value \\J~q 1 - A simple analysis 
shows that the mean quark energy (E q ) is related to the mean gluon energy (E g ) calculated 
in ||. Although this result is quite plausible at first sight it is nevertheless surprising since 
the calculation of the gluon's mean energy (E g ) involves tree graph contributions only 
whereas one needs to bring in loop contributions for the determination of the quark's 
mean energy (E q ). In Sec. 3 we check on the rigidity of the back-to-back top quark pair 
configuration against gluon radiation by looking at opening angle distributions and their 
mean values. We determine the differential cos ^-distribution where #12 is the opening 
angle of the top quark pair and compare to the corresponding distribution in the bottom 
quark case. We then go on to calculate the mean deviations of cos 9 12 and sin #12 from 
the anticollinearity limits cos 9u = — 1 and sin 9 12 = 0. In Sec. 4 we study mean values of 
the longitudinal energy El and the transverse energy E? of the top quark. The leading 
contributions to these measures are given by the mean energy (E q ) calculated in Sec. 2 and 
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by the mean opening angle measures (1 + cos 6*12) and (sin #12) calculated in Sec. 4. We also 
numerically evaluate the mean transverse momentum (py) of the top quark relative to the 
antitop direction and the mean of the opening angle itself. Sec. 5 contains our summary 
and conclusions. In Appendix A we list explicit formulas for the full 0(a s ) heavy quark 
production cross section which enter in the evaluation of the mean values as denominator 
functions. In Appendix B we list a number of basic rate functions which appear in our 
analytical results. In Appendix C we present analytical results on first moment integrals 
that are needed for the evaluation of the mean values of (E q ) and (1 + cos #12) for scalar 
and pseudoscalar current- induced top pair production. 



g 11 H 1 (y,z)+g 12 H 2 (y,z)) (1) 



2 Quark energy spectrum and (E q ) 

To begin with we express the differential heavy quark production cross section in e + e~ — > 
q(Pi)<l(P2)g(P3) with q = pi + p 2 + P3 in terms of the relevant components of the hadron 
tensor H< = {-g^ + q^q u /q^Hj^ 

da a 2 
dy dz ASixq 2 

where y — 1 — 1p\qjq 2 and z = 1 — 1p 2 qjq 2 are fractional energy loss variables for the 
heavy quark and antiquark. We shall also use the abbreviation £ = 4m 2 , /q 2 . Another 
useful variable is the Born term velocity of the quark given by v = y/l — £. We closely 
follow the notation employed in |J. Thus the superfix i — 1,2 stands for the two current 
combinations H 1 = \(H VV + H AA ) and H 2 = \(H VV - H AA ). In Eq. (§ we have already 
averaged over the relative beam-event orientation. In the terminology of || the full hadron 
tensor components would thus carry the additional index (U + L) relevant to the total rate 
(U: unpolarized transverse, L: longitudinal). However, since we are only concerned with 
the total rate in this paper we have consistently dropped this index. The parameters gn 
and #12 specify the electro-weak model dependence and are given by 

9n = Q) ~ 2Q f v e v f Rexz + (v 2 e + a 2 e )(v 2 f + a 2 f )\xz\ 2 , 

912 = Q 2 f -2Q f v e v f Re X z + (v 2 e +a 2 e )(v 2 f -a 2 f )\xz\ 2 , (2) 

where, in the Standard Model, xz(q 2 ) = gM\q 2 {q 2 - M§ + iMzTz)- 1 , with M z and 
T z the mass and width of the Z° and g = G F (8v / lvra)- 1 w 4.49 • 1(T 5 GeV' 2 . Q f are 
the charges of the final state quarks; v e and a e , Vf and a/ are the electro- weak vector 
and axial vector coupling constants. For example, in the Weinberg-Salam model, one has 
v e = — 1 + 4 sin 2 9 W , a e = — 1 for leptons, Vf — 1 — | sin 2 9w, a/ = 1 for up-type quarks 
(Qf = |), and Vf — —1 + Isin 2 ^, cbf = —1 for down-type quarks (Qf = — |). In this 
paper we use Standard Model couplings with sin 2 9w = 0.226. 

The hadron tensor components H l can be obtained from the relevant tree-level Feynman 
diagrams by the above projection H l = (—5^ + q^q v /q 2 )H. \ v . They are given by 



H\y,z) = N 
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where N = 4tcq: s NcCf- The hadron tensor component H 2 (y,z) features the typical chi- 
rality factor £ which guarantees that H vv = H AA in the limit v — > 1 or m q — > 0. 

In order to obtain the energy spectrum of the quark we integrate over the fractional 
energy loss variable z of the antiquark. The integration limits are given by 



z± 
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The integration is straightforward and results in the hadron tensor components 
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In Fig. la we show the resulting ^/-spectra for bottom and top pair production at ^fq 1 = 
500 GeV using m& = 4.83 GeV, m t = 175 GeV and a running a s with a s (m 2 z ) = 0.1175 
(rif = 5 below tt-threshold and n/ = 6 above tt-threshold). Both spectra are strongly 
peaked towards low values due to the infrared (IR) singularity at y — 0. The differential 
rate for bottom pair production is much higher than that of top pair production because 
the phase space is much larger for bottom pair production. 

Next we determine the mean energy of the quark. Up to 0(a s ) the mean energy of the 
quark is defined by 



(E q ) 



— \J q 2 a (Born) 



da(tree) 1 r^r 

E q — - — - — dydz + -yq z cr(loop) 
dydz I s 



a(Born) + a(a s 



(8) 



where a(a s ) is the 0(a s ) cross section a(a s ) = cr(tree) + a (loop). Note that we define our 
mean energy relative to the full 0(a s ) cross section and not to the Born term cross section 
as done in ||. We shall come back to this point later on. In order to make the paper self- 
contained we list the Born term and 0(a s ) cross sections appearing in the denominator of 
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Eq. (H) in Appendix A. We have checked that the infrared (IR) singularity resulting from 
the tree-graph integration in the numerator of Eq. (§) cancels against the IR- singularity of 
the loop contribution, i.e., including the loop contribution, the mean energy of the quark 
is a manifestly IR-safe measure. In fact Eq. @ can be rewritten into a more compact form 
by using the fractional energy loss variable y with E q = — V)- One has 



f dcr(tree) 

/ V— j— 5 — dydz 
J dydz 



(y) = (1 - 2£ g / V<f> = , H M 7 r - (9) 

The calculation of the mean (y) is simpler since one only needs to compute the tree- 
graph contribution. This is plausible since by taking the first moment of the differential 
tree-graph cross section with regard to y the IR-singular piece in the tree-graph cross 
section is cancelled. In the following we shall refer to the mean of variables that vanish at 
the IR-point as tree-graph IR-safe measures. For tree-graph IR-safe measures such as (y) 
the perturbation series starts at 0(a s ). At this order the inclusion of the full 0(a s ) total 
cross section in the denominator of Eq. (H) is no longer mandatory. By keeping the full 
0(a s ) contribution in the denominator one partially includes higher order effects in the 
evaluation of (y). In order to be definite, we shall always present our results on tree-graph 
IR-safe measures normalized to the full 0(a s ) cross section. 

We give the results of the first moment integration of the numerator of Eq. (|^) in terms 
of the first moments of the two relevant hadron tensor components defined by 

H*M = JyH\y,z)dydz. (10) 

We obtain 
H i[y] = N 



-"-(256 - 224£ + 120£ 2 - 8£ 3 - 5£ 4 )t 3 - 7^(704 + 336£ + 34£ 2 + 15£> 
24 36 



and 

#%] = —(128 - 120£ + 36£ 2 + 5£ 3 )t 3 - —(280 - 118f - 15£ 2 )v 
24 36 
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(12) 



Our results for the y-moments are simply related to the fractional gluon energy x- 
moments H VV W = H^+H 2 ^ and H AA ^ = H^-H 2 ^ calculated in § (x = 2p 3 q/q 2 = 
2E g / \fq I ) as can be seen by the following reasoning. Due to CP-invariance of the under- 
lying interaction, the y-moments for the quarks and the z-moments H 1 ^ for the 
antiquarks are identical, i.e. H 1 ^ = H^ Z K From the kinematical relation x = y + z one 
then obtains 

H i[y] = _ H i[x) ( 13 ) 

as already noted in ||. Our moment results agree with the results given in ||. 

In Fig. 2a we plot the mean value of the fractional energy loss variable (y) = (1 — 
2E q / ^/q 2 ) for bottom and top pair production as a function of y/q 2 . We first discuss the 
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bottom pair production case in Fig. 2a where the mean value (y) shows a sharp rise from 
threshold and then levels off at higher energies. To understand the high energy behaviour 
we consider the limiting behaviour of the y- moments H 1 ^ as v — > 1. In the high energy 
limit v — > 1 or, equivalently, in the limit m q —>■ we obtain 
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H 2[y] -> 0. (14) 



One notes that if 1 ^] becomes mass singular in this limit due to the presence of the 
quark/ ant iquark-gluon collinear singularity. There is no inconsistency in this singular 
behaviour from the physics point of view because the measure (y) is no longer defined in 
this limit. The reason is that the energy of the quark cannot be separately measured in the 
collinear mass zero case. In the corresponding collinear massive configuration the energy of 
the quark can in fact be separately determined by e.g. a time-of-flight separation. In order 
to define a physically meaningful (y) in the m q — > limit one has to introduce an angular 
cut-off. For very small values of the angular cut-off one would have to exponentiate the 
resulting large logarithmic factor taking into account multiple gluon emission. Returning 
to Eq. (14) the vanishing of H 2 ^ is expected since H vv — > H AA for v — > 1 or m q — > as 
mentioned before. The limiting behaviour of (y) following from Eq. (|L1]) can be seen to be 
given by 

g.321n(4/fl-176/3 - 

7T 36(1 + « s /7r) 

The logarithmic growth in Eq. ([15]) resulting from the mass singularity is off-set by the 
logarithmic fall-off of a s . Looking at Fig. 2a the mean value (y) is still slightly rising at 
\[q 1 = 1000 GeV showing that the approach of (y) to its asymptotic value is not very fast. 

Turning to the top pair production case in Fig. 2a we first determine the threshold 
behaviour of (y). In the nonrelativistic limit v — > we obtain 

H m ^. Nv b + (v 7 ), H 2 ^ - —Nv 5 + 0(v 7 ). (16) 
5 5 

The limiting behaviour of the H 1 ^ can be understood in the following way. At threshold the 
cross section factorizes into the Born term cross section e + e~ — > qq and a gluon production 
piece. For the Born term cross section the s-wave vector current contribution dominates 
over the p-wave axial vector current contribution leading to H 1 ^ ~ H 2 ^ close to threshold. 
The leading v 5 -behaviour of the s-wave vector current contribution has its origin in the 
following factors: there is one power of v from the flux factor and four powers of v from 
integrating the y- weight over gluon phase space. We do not dwell on the nonleading 0(v 7 )- 
contributions in Eq. ( |i~6"D but mention that there is some interesting physics concealed in 
the nonleading 0(v 7 )-contributions as discussed in 0. Finally, when calculating (y) one 
divides by the Born term cross section resulting in an overall w 4 -threshold behaviour. From 
what has been said about the low and high energy behaviour of (y) it is clear that close 
to tt-threshold the vector current contribution is the dominant contribution whereas at 
high energies the vector and the axial vector currents contribute equally to (y) . The same 
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statement is true for the other linear mean values calculated in the following sections. We 
want to mention that the high energy limit (v — > 1) and nonrelativistic (v — > 0) limits of 
(y) presented in this paper are in complete agreement with the results given in |J. 

The power behaved dependence of (y) on y/q 1 close to top-antitop threshold is clearly 
visible in Fig. 2a. Away from threshold (y) starts turning over at around 500 GeV in 
its slow approach to asymptotia. At 500 GeV the fractional energy loss of the top quark 
amounts to less than 1% while at 1000 GeV the fractional energy loss is 3.77% (see Table |l|). 



mean value 


yfq* = 500 GeV 


Vq 1 = 1000 GeV 


(y) 


0.71% 


3.77% 


(1 +COS#i 2 ) 


0.41% 


1.99% 


(y(l + cos0i 2 )) 


0.06% 


0.57% 


(y L ) 


1.06% 


5.19% 


(sin 12 ) 


2.07% 


7.30% 


(ysin 12 ) 


0.19% 


1.23% 


(Vt) 


1.88% 


6.06% 


(x) 


1.43% 


7.53% 




1.18% 


5.46% 




1.25° 


4.62° 




0.490 


0.123 


V 


0.714 


0.997 



Table 1: Mean values of different kinematical quantities in e + e — > ttg for yfq 1 = 500 GeV 
and 1000 GeV. We also list the values of £ and v = y/1 — f. 



3 Opening angle distribution 

and the mean values (1 + cos On) and (sin #12) 

Next we determine the differential distribution in the cosine of the opening angle cos6 , 12 
and the mean deviation from the anticollinearity limit cos6 l 12 = — 1. The central relation 
is the relation between cos Q\i and the (y, z)-variables which reads 

yz + y + z-l + £ 
cos 6>i2 = , — , (17) 

;i-y) 2 -W(i-*) 2 -£ 



For later purposes we need to solve Eq. (^7|) in terms of the y-variable. One obtains a 
quadratic equation in y which can be solved to give 

(1 - z)(l - COS0 12 ) + z((2 - z) cos 2 gig - (z + Q) 
y± ~ (l + ^) 2 -cos 2 12 ((l-^) 2 -O 



cos0W(l - zf - UAz* - £((1 - zf - 0(1 - cos. 



12; 



(l + z) 2 -cos 2 ^ 12 ((l-^) 2 -0 * v ' 

First we discuss the differential cos ^-distribution. With the help of relation ( jETf) the 
differential (y, ^-distribution given in Eq. (|IJ) is transformed into 

da da dy 



d cos 6 12 dz dy dz d cos 12 
where the partial derivative is given by 

dy = ((i-^) 2 -0 1/2 ((i-y) 2 -0 3/2 f2m 

«9cosfl 12 (i- y )(2z + £)-Z(z+l) ' 1 J 

After substituting for y one can then integrate Eq. ([B5]) with regard to z. Bearing in 
mind that the two solutions of the quadratic equation Eq. ([TBI) y± = y±(z, cos 6*12) have to 
be substituted in the hadron tensor components H l (y,z) and the Jacobian Eq. fl20| ) it is 
evident that it is a cumbersome task to do the requisite z-integration analytically. Instead 
the integration will be done numerically. 

A closer look at phase space is of great help in disentangling the ^-integration limits 



and in determining which of the two solutions y± = y±(cos8i 2 , z) of Eq. (18) have to be 
substituted for the y- variable in Eqs. (|3]), (|4]) and fl20"|). In Fig. 3 we show a plot of the (y, z)- 
phase space for m t = 175 GeV and a/*? 2 " = 500 GeV. Apart from the phase space boundaries 
defined by cos 6*12 = ±1 we have included the contour lines cos 6*12 = 0, ±0.7. All contour 
lines of constant cos#i 2 must intersect at the points (ija, za) = (1— v^l - \/£)/(2— Vf)) 
and (yB, zb) = (Vf(l ~~ \/£)/(2 — 1 ~~ V£) where the antiquark and the quark resp. are 
at rest. The reason is that at these two points the opening angle 8yi is no longer defined. 
It is evident that the upper limit of the z- integration is always given by zb = 1 — v^- 
The lower z-limit depends on whether cos 812 is positive or negative. For positive values of 
cos 812 the lower limit is given by za = VX(l — v^D /(2 — vX) • F° r negative values of cos 6*12 
the boundary curve zc = 2 ; c(cos6 l i 2 ) can be determined by first solving for the minima 
of the z = z(cos8i2,y) contours in the (z, ?/)-plane and then by inserting the y-loci of the 
minima into Eq. (|TTD. For the position of the z- minima one finds 

V= Z 4^- (21) 
2z + y K ' 



After substituting (|2TJ) into Eq. ( |17D one obtains a quartic equation for zc = zc(cos#i 2 ) 



From the four solutions of the quartic equation the one relevant to the case at hand is 
given by 



-£(1 - cosM + J£(l - cos0 12 )(4(l - + e 2 (l " cos# 12 )) 
4 — 4(1 — cos6>i 2 ) 

In Fig. 4 we show a plot of the phase space in the (z, cos #i 2 )-plane with boundary curves as 
determined by the above discussion. We have also indicated which of the two solutions of 
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the quadratic equation y± = y(z, cos #12) of Eq. ([Tj|) have to be used in the two respective 
phase space regions 1 < cos #12 < and < cos 6*12 < — 1. 

In Fig. lb we finally show the differential cos ^-distribution for bottom and top pair 
production again at the center of mass energy of 500 GeV. As expected the distributions 
are strongly peaked towards the IR-point at cos #12 = — 1. The differential bottom cross 
section dominates over the differential top cross section because of its larger phase space. 

Next we determine the mean deviation from the anticollinearity limit cos #12 = — 1. To 
0(a s ) the mean of 0086*12 is defined by 



(cos 6»i 2 ) = 



—a(Born) + [ cos9 12 — - — : — -dy dz — a (loop) 
J dydz 



a(Born) + cri 



a.s 



(23) 



Again we have checked by explicit calculation that the IR-singularity from the tree graph 
integration is cancelled by the corresponding IR-singularity of the loop contribution. Sim- 
ilar to Eq. (|9p one may banish the IR-singularity in the tree graph contribution by taking 
the moment with regard to the weight factor (1 + cos 612). In explicit form one adds and 
subtracts unity in the tree-graph integral, cos 9 12 = 1 + cos 9 12 — 1. One can then rewrite 
Eq. (|D as 

da (tree) 



(1 + cos #12) 



/ (1 + cos 0i 2 ) , , dydz 
J dy dz 



(24) 



a(Born) + a(a s ) 

What is needed are the first moment integrals H 1 ^ taken with regard to the moment 
variable t = l+cos^ for the numerator functions. Since (l+cos^) is a tree-graph IR-safe 
measure on needs to consider only the tree-graph components of the 0(a s ) hadron tensor 
as explicated after Eq. @. The requisite moment integrations can be done analytically. In 
fact similar integrations appear in the calculation of longitudinal spin-spin correlations || 
H , and some of the analytical results can be taken from these papers. One obtains 



N 



16^ - 62£ + 56^ + 3£ 2 + 3£ 2 ^ + 4(4 - 5£ - 50)- 

V V v v 



+ (128 - 128£ - 120 + 390 - 30 - 32(4 - £)v 3 ) ^ 



+2(4 - 0(2 - 0(*9 - h*) - (64 - m - 116£ 2 + 31f 3 - 3£ 



2 , JI( , _ 

2v 2 



-4(4 + 0(1 - 20*14 + 2(16 - 40£ + 310 - 10f 3 ) 



#2[i] 



v 

98 - 68 x /^ - 3£ - 30/^ - 4(7 - £)- + 3(32 - 250 + £ 3 - 96w 3 )-% 
v v v 4^2 



(25) 



+6(2 - 0(t 9 - tie) + (208 - 168£ + 350 - 3£ 3 ^ 13 



-4(13 + 30^14 + 6(4 - 13f + 80)% 



2v 2 



(26) 
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The rate functions U (for i — 3, 9, 13, 14, 15, 16) appearing in Eqs. ( |25| ) and ( p6|) are listed 
in Appendix B. As will be discussed in Sec. 4, (1 + cos6 l 12 ) together with (y) determine 
the leading contribution to the mean of the longitudinal energy (El). 

In Fig. 2b we show a plot of the v^-dependence of (1 + cos 6*12) for bottom and 
top production. In the bottom case (1 + cos #12) shows a quick rise from threshold and 
very quickly reaches its asymptotic value away from threshold. The asymptotic value is 
determined by the limiting behaviour of the t-moments in the high energy limit v — > 1. 
These are given by 

H 1[t] -> 24iV, H 2[t] -> 0, (27) 

which translates into 

(l + cosfl 12 >^ — ; , 2 / • (28) 

In Fig. 2b we show a plot of the energy dependence of (1 + cos #12) for bottom pair pro- 
duction. After a steep rise from threshold the mean value (1 + cos^) starts to decrease 
again at around 200 GeV. The logarithmic fall-off which sets in beyond 200 GeV is due 
to the combined effect of (1 + cos 812) reaching its asymptotic value and the logarithmic 
running of the coupling constant a s . 

Turning to top pair production in Fig. 2b we first discuss the threshold behaviour of 
(1 + cos 812). In the nonrelativistic limit v — > the threshold dependence of the moments 
fl*m i s given by 

H m ^Nv 5 + 0(v 7 ), H 2 ^ - — Nv 5 + 0(v 7 ). (29) 
75 v ' 75 K ' K ' 

One thus again has a w 4 -threshold behaviour for the mean (1 + cos 612). The t> 4 -threshold 
behaviour can clearly be discerned in Fig. 2b. At around 600 GeV (1 + cos 812) starts 
turning over in its approach to its asymptotic value. Numerical values of (1 + cos #12) for 
top pair production are given in Table [I]. 

Next we determine the mean value of sin #12 which, according to the discussion in Sec. 4, 
gives the leading contribution to the mean transverse energy of the quark. The mean value 
(sin 812) provides another measure of how much gluon radiation distorts the lowest order 
back-to-back configuration of the quark-antiquark pair. First we express sin# 12 in terms 
of the (y.z)-phase space variables. One has 



. jAyzil - y - z) - t{y + zf 
sin 012 = . ; (30) 

:i-y) 2 -W(i-^) 2 -^ 



Similar to (1 + cos 812) the mean value of sin# 12 is determined in terms of the 0(a s ) tree 
graph contribution alone since sin# 12 vanishes at the IR-singular point. In the terminology 
of the previous sections (sin 812) is a tree-graph IR-safe measure. The necessary integrations 
have been done numerically. The top production results are shown in Fig. 5b where we plot 
(sin #12) as a function of the center of mass energy \f(p. A comparison with Fig. 5a shows 
that (sin #12) is smaller than (1 + cos #12) over the whole range from top-antitop threshold 
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to 1000 GeV. Qualitatively this is not difficult to understand since close to the IR-point 
cos 9 12 = — 1 where the cross section is largest the weight function sin^ 12 = \/l — cos 2 6>i 2 
clearly dominates over the weight cos6 l 12 when considered as functions of cos #12. In fact 
sin #12 has an infinite slope at the IR-point compared to the unit slope of cos6 , 12 . One 
therefore expects (1 + cos #12) < (sin #12) which is quantitatively borne out as Figs. 5a 
and 5b show. Numerical values for (sin6 l i 2 ) at the two center of mass energies 500 and 
1000 GeV are listed in Table @. 

Of interest is also the mean value of the opening angle 9u or its complement, the 
acollinearity angle 9\i = 180° — 9rz- In Fig. 6 we show a plot of the energy dependence 
of the mean of the acollinearity angle (9u)- (9 12) rises smoothly from threshold to a 
value of 1.25° at 500 GeV and further to 4.62° at 1000 GeV. A good estimate of (0 12 ) 
can be obtained by inverting the corresponding numerical values for (sin #12) in Table [I] 
which gives arcsin((sin #12)) = 1-19° and 4.19°, respectively. Using instead (l + cos#i2) for 
the corresponding estimate one obtains 5.16° and 11.45° which considerably overestimates 
the true values of (#12)- The reason that (sin6 l 12 ) gives a good estimate of the mean of 
the acollinearity angle (#12) can again be traced to the fact that the weight factor sin 9 12 
(considered as function of cos 6*12) has a strong support close to the IR-point cos 9 12 — — 1. 



4 Mean values of the transverse 
and longitudinal energy 

Up to now we have been separately considering the energy loss of the quark and the 
deviation of the acollinearity angle 9 12 from the anticollinearity limit 9 12 = 180° due to 
gluon radiation. From the physics point of view it is also interesting to consider the 
energy loss of the quark projected in the longitudinal direction as well as the build-up 
of transverse energy resulting from gluon radiation. Here the longitudinal and transverse 
energies are defined relative to the momentum direction of the antiquark and are thus 
defined by El = —E q cos9i 2 and E? = E q sm9i 2 . Again we shall calculate mean values 
of these two quantities as a measure of the importance of gluon radiation. As before we 
rewrite (El) and (Ej) in terms of tree-graph IR-safe measures. For the mean longitudinal 
energy one has 



(e l ) = \>J? 



(1 + cosfl 12 > - (y) + (y(l + cos9 12 )) 



(31) 



or, written in terms of the fractional longitudinal energy loss variable yL = 1 — 2EL/yfq I , 
one has 

(y L ) = [(1 + cos 9 12 ) + (y) -{y(l + cos 12 ))] . (32) 

The mean values appearing in Eq. ( |3"2"D must satisfy three geometrical inequalities which 
follow from the geometrical inequalities y < yL, (1+cos 6 I 12 ) < yL and y—y( 1+ cos 9 12 ) < yL- 
Since these inequalities are true on an event-by-event basis they must also hold for their 
means. One can check that the relevant numerical entries in Table [l| satisfy the three 
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inequalities. On top of the geometrical inequalities one may classify the mean values 
of the different weight variables in Eq. ([32]) in terms of powers of smallness since the 
differential cross section is strongly peaked towards the IR-region. Thus one expects that 
the leading contributions to (El) come from the linear terms (1 + cos6 l i2) and (y) whereas 
the quadratically small term (y(l + cos #12)) is expected to be a nonleading effect. This 
qualitative picture is borne out by the numerical values in Table p] and by the curves 
in Fig. 5a where the energy dependence of the above mean values are combined in one 
common plot. 

The mean of the fractional transverse energy y? = 2Et/\/~(P can be written as 

(Vt) = [(sin6»i 2 ) - (ysme 12 )\. (33) 

One now has the geometrical inequality (y T ) < (sin 12 ) which can be seen to be satisfied 
by the relevant numerical entries in Table |l|. The leading contribution is expected to arise 
from the linear term (sin 0i 2 ) whereas the contribution of the term (ysinO^) is nonleading. 
The numerical values in Table |] as well as the relevant curves in Fig. 5b confirm this 
picture. 

Even though we do not present numerical results for the bottom production case we list 
the high energy limits of the mean values (y(l + cos0i 2 )), (sin0i 2 ) and (y sin 6*12) appearing 
in Eqs. ( |3"2"D and ( ]3lf ) for the sake of completenes. Including also the result Eq. flZgp one 
has 

(1 + cos 12 ) - 
(2/(1 + cos 0i 2 )) - 
(sin 0i2 ) - 



a s 


2 


71 


1 + a s /V 


a s 


20 


71 


27(1 + a s /n) 


a s 


IOtt 


71 


9(1 + « s /tt)' 


a s 


71 


7i 


2(l + a s /7r)' 



(ysinM - ~ on , /v (34) 
7r 2(l + a s /7r) 

Where applicable the limiting values can be seen to satisfy the geometrical inequalities 
written down before. The limiting values show the same hierarchy in terms of powers of 
smallness as discussed before for top production below 1000 GeV. However, the hierarchy 
is not as pronounced as in the top production case below 1000 GeV. As concerns the 
limiting values in Eqs. (0) we want to add two remarks. In the mass zero limit treated 
above the transverse energy of the quark is equal to its transverse momentum. Also the 
mean of the longitudinal energy (El) and the longitudinal fractional energy loss (yi) in 
Eqs. (|3l|) and ([32|) are mass singular as manifested by the contribution of (y) to (El) and 
(yL) which by itself is mass singular (see Sec. 2). 

As Table shows, gluon radiation causes an average energy loss of 1.06% of the top 
quark's energy in the longitudinal direction and an average energy gain of 1.88% in the 
transverse direction at this energy. At 1000 GeV the corresponding percentage figures 
are 5.19% and 6.06%. This has to be compared with the total energy loss of 0.71% at 



12 



500 GeV and 3.77% at 1000 GeV as calculated in Sec. 2. The gain in transverse energy is 
larger than the loss of longitudinal energy in this energy range where this effect becomes 
smaller as the center of mass energy increases. Qualitatively this behaviour may be un- 
derstood from what is referred to as the "dead-cone" effect of gluon emission from a heavy 
quark fTOf . Gluon emission is suppressed in an angular cone of size 9 ~ = 2m q /\f(p 
around the quark-ant iquark direction, i.e. the emitted gluon has a strong momentum com- 
ponent transverse to the quark-antiquark direction. Because of momentum conservation 
the transverse momentum of the gluon has to be balanced by the transverse momentum 
of the quark. A good estimate of the fractional transverse momentum of the heavy quark 
in the energy region under consideration can be obtained by rescaling yx by the factor 
v = p/E = 2p/y/q I where p is the lowest order momentum of the quark (see Table |l|). 
Taking the rescaled values of (i/t) and comparing it to the mean fractional total energy of 
the gluon (x) in Table |I]one concludes from momentum conservation that, on average, the 
gluon must have a large transverse momentum or energy component. By comparing the 
relevant numbers at 500 and 1000 GeV one notes that this effect becomes smaller as the 
energy is raised in agreement with the shrinking of the "dead-cone" predicted in []TU| . 

In order to get a more quantitative handle on the "dead-cone" effect we show a plot of 
the g 2 -dependence of the mean fractional gluon energy (or momentum) (2E g / y/(p) = (x) 
and the mean fractional transverse momentum (2px/ \fq I ) of the top in Fig. 7. In the whole 
range from threshold up to 1000 GeV the mean transverse momentum of the top quark is 
almost balanced by the mean total momentum of the gluon where the approximate equality 
of the two means degrades as the energy increases. This proofs that the emitted gluon 
has a dominant transverse momentum component, and that the preference for transverse 
emission is degraded as the energy increases, just as expected from the "dead-cone" effect. 



5 Summary and conclusions 

We have discussed in some detail how gluon radiation affects the lowest order back-to- 
back configuration of pair produced top quarks in e + e~ annihilation. We have introduced 
various measures that serve to quantify gluon radiation effects in this reaction and we have 
studied their energy dependence. These measures are related to the overall energy loss and 
the energy loss in the longitudinal direction, or, in the case of the transverse energy, to 
the energy gain in the transverse direction. We have also introduced two angular measures 
related to the opening angle #12 between quark and antiquark. 

When possible we have given analytical expressions for the energy and momentum mea- 
sures (E q ), (E T ), (Ei) and (px) and for the angular measures (1 + cos #12) and (sin^ 12 ). 
We have provided numerical results for all these measures in the energy range from close 
to top-antitop threshold up to 1000 GeV. For some of the measures we have compared 
the top-antitop results to the corresponding figures for bottom pair production including 
a discussion of the high energy limit v — > 1 limits of the various measures. The numerical 
predictions and their energy dependence can directly be compared to experimental data 
on top pair production at the proposed future linear colliders. The values of these mea- 
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sures also serve to quantify the effect of gluon radiation on the lowest order back-to-back 
configuration of the produced top quark pairs. 

Compared to the lowest order energy of the top E q = \\fq I the mean loss of longitudinal 
energy is smaller than the mean gain in transverse energy in the range of cm. energies 
from threshold to 1000 GeV . At 500 GeV the mean transverse energy exceeds the mean 
longitudinal energy loss by about a factor of two. In this sense the top quark is more rigid in 
the longitudinal direction than in the transverse direction. This behaviour can be traced 
to an effect which is referred to as the "dead-cone" effect in the literature. In massive 
quark pair production with additional gluon radiation the gluon tends to be emitted away 
from the heavy quark and antiquark directions. Since the mean momentum of the gluon is 
comparatively large ((x) = 2(y)) and since the transverse momentum of the gluon has to be 
balanced by the transverse motion of the quark (or antiquark) one obtains comparatively 
large values for the mean transverse energy and momentum of the quark (or antiquark). 



Appendix A: 0(a s ) cross section 



The complete 0(a s ) cross section for heavy quark production in e + e~ — > qqg has been 
calculated before in || |7|]. In order to make the paper self-contained we list the results in 
this Appendix. For the Born term cross section a (Born) one obtains 



where 



a 1 (Born 



a(Born) = [gua 1 (Born) + gi 2 cr 2 (Born) 
na 2 vNc 



-(4-0, a 1 (Born) 



7ra 2 vNc 



3£- 



(Al) 



(A2) 



3q 2 ^ 3q 2 

The 0(a s ) tree- and loop-contributions are given in terms of the hadron tensor components 
H % defined in Eq. ([TJ). One has 

|(4-o(2-o«+^i 

-2(4 - 0(2 - 0(*8 - t 9 ) - 4(4 - O^io + 2t 12 )l , (A3) 



H l (a s ) 



N 



\(A - 0(2 - i)v + -(192 - 104£ - A? + 3£ 3 )*3 



H 2 (a K 



;i8_£) v + 7 (24-8£-r)t3 

2*i 2 ) 

The rate functions *j (for i — 3, 8, 9, 10, 12) are given in Appendix B. 



2 y " / 4 
6(2-0(*8-* 9 )-12^(*i 
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Appendix B: Rate 

The analytical expressions of 
involve a number of basic rate 



integrals and first 

the 0(a s ) cross section and 
functions which are collected 



moment integrals 

the mean value (1 + 0086*12) 
in this Appendix. 
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Appendix C: Scalar and pseudoscalar currents 

In this Appendix we present expressions for the hadron tensor components and their first 
moments which allows one to calculate the mean values (y) and (1 + cos# 12 ) for the scalar 



KuSS 



H ™ - 



(S) and pseudoscalar (P) current contributions H 1 = ^(H ss + H pp ) and H 2 
H pp ). These expressions would be needed for the assessment of gluon emission effects in 
the heavy quark production process e + e~ — > qq(g) mediated by scalar (S) or pseudoscalar 
(P) particles as resulting from e.g. Higgs exchange. One has 
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Figure Captions 



Fig. 1: Differential (a) y- and (b) cos # 12 -distribution for e + e — > bbg and e + e — > 
at = 500 GeV (m* = 175 GeV, m b = 4.83 GeV) 

Fig. 2: Mean values of (a) the fractional energy loss variable y of the bottom and 
top quark and (b) the variable t = 1 + cos 9 12 as functions of sfq 1 for bottom 
and top quark pair production (m t = 175 GeV, nib = 4.83 GeV) 

Fig. 3: (y, z) -phase space for e + e _ — > tt<? (m t = 175 GeV, = 500 GeV). Also 
shown are contour lines for constant values of cos # 12 . To the right and to 
the left of the zero contour line cos#i 2 = are the phase space regions 
< cos 9 12 < 1 and —1 < cos#i 2 < 0. All contour lines of constant cos#i 2 
intersect at points A and B. Point C denotes the minimum of the cos #12 
contour lines for negative values of cos#i 2 . The IR-point is at the origin. 

Fig. 4: (z, cos#i 2 )-phase space for e + e~ — > tig (m t = 175 GeV, y/q 1 = 500 GeV). 

The IR-point is at the lower left corner. The dark-shaded y + - and the light- 
shaded y_-regions correspond to the two solutions of the quadratic equation 
V± = y±(z, cos 6>i2 ) 

Fig. 5: Comparison of the mean values of (a) y^ — 1 — lEi^j^fq 1 , y, 1 + cos6 l 12 , and 
2/(1 + cos #12) and (b) yx = 2_& r /y / g 2 ', sin# 12 , y, and ?/sin# 12 for top pair 
production in e + e~ — > tig as function of y/q 1 

Fig. 6: Mean value of the acollinearity angle 9± 2 in e + e~ — > tig as function of \fq I 

Fig. 7: Mean values of the fractional gluon energy x and fractional transverse mo- 
mentum 2pT/yfq I of the top quark in e + e~ — > ttg as functions of 
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